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I. INTRODUCTION 

Test particle motion in stochastic velocity fields is a 
generic problem in various topics of fluid and plasma 
turbulence or solid state physics [1]- [4]. This problem 
was very much studied and several analytical approaches 
were developed. Most of them are based on Corrsin [5], 
[1] and direct-interaction [6], [7] approximations. 

In this context, particle motion in 2-dimensional 
divergence-free velocity fields represents a special case. 
Kraichnan has shown for the first time, in a study based 
on numerical simulations [8], that the existing analyti- 
cal methods are not adequate for this type of problems. 
The cause of this anomaly is the trapping of the parti- 
cles, which appears in such velocity fields when they have 
slow time variation. The trapping consists in trajectory 
winding on almost closed paths. A typical trajectory has 
a complicated shape with such localized trapping events 
separated by long jumps. Consequently, the probability 
distribution function is non-Gaussian. Thus Corrsin and 
direct-interaction approximation which are based on this 
hypothesis are not adequate for this specific case. A more 
recent analysis of the effects of trapping is presented in 
Ref. [9] where a non-Gaussian peaked distribution of the 
displacements and a long negative tail in the Lagrangian 
velocity correlations are evidenced for numerically calcu- 
lated trajectories. 

This special problem of diffusion in 2-dimensional 
divergence-free velocity fields describes for instance the 
transport in turbulent magnetized plasmas or in incom- 
pressible fluids. It was studied especially by means of di- 
rect numerical simulations ( [10] and the reference there 
in) or on the basis of simplified models [11], [12]. There is 
also a qualitative theoretical estimation of the scale law 
for the asymptotic diffusion coefficient [13] based on an 
analogy with the percolation process in stochastic land- 
scapes. The case of collisional particle motion in such 
static velocity fields was analyzed by means of the renor- 
malization group techniques ( [2] and the references there 
in) and the asymptotic time behavior of the mean square 
displacement was determined. The evolution of the diffu- 



sion process is determined only in [14] where a statistical 
approach, the decorrelation trajectory method, is devel- 
oped. It yields analytical expressions for the time depen- 
dent diffusion coefficient D{t) and for the correlation of 
the Lagrangian velocity L{t), that are qualitatively valid 
for the whole range of the Kubo number (see next Section 
for the definitions). The basic idea consists of determin- 
ing the Lagrangian velocity correlation by means of a set 
of average Lagrangian velocities estimated in subensem- 
bles of realizations of the stochastic field. This method 
could be extended to more complicated physical systems 
which contain particle collisions [15], average velocities 
[16] or a supplementary component of the motion perpen- 
dicular to the 2-dimensional plane [17], [18]. It was shown 
that the presence of trapping strongly influences the dif- 
fusion coefficients and their scaling laws determining a 
rich class of anomalous diffusion regimes. These stud- 
ies have shown that the decorrelation trajectory method 
provides a qualitatively good description of the trapping 
process. However, due to the rather strong approxima- 
tion introduced in this method (see Section III), there 
are several qualitative aspects that are not well described 
[19], [20]. They are related to trajectory fluctuations and 
their correlation with the stochastic velocity. 

The above results concern the effect of trapping on 
the individual trajectories. The trapping has also col- 
lective effects. It determines coherence in the stochastic 
motion in the sense that bundles of neighboring trajec- 
tories form localized structures similar with fluid vor- 
tices. The formation of these structures strongly influ- 
ences the transport. The aim of this paper is to study 
the statistical characteristics of these trajectory struc- 
tures determined by the intrinsic trapping appearing in 
2-dimensional divergence-free stochastic velocity fields. 
The average, the dispersion and the probability distribu- 
tion function for the trajectories in such structures are 
determined. The evolution of these statistical quanti- 
ties saturates showing that they do not contribute to 
the asymptotic transport. The statistical evolution of 
the distance between two neighboring trajectories is also 
studied and its mean, dispersion and probability distri- 
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bution arc determined. We show that for the trapped 
trajectories in a static velocity field, the average decays 
to zero and the dispersion remains very small for long 
time and eventually saturates. A very strong clump ef- 
fect is found for the trapped trajectories. This demon- 
strates the existence of trajectory structures and their 
effect of strongly reducing the relative transport. The 
later is produced only by a small part of the stochas- 
tic trajectories which are not contained in these vortical 
trajectory structures. 

The method developed for this study is a semi- 
analytical statistical approach which extends and im- 
proves the decorrelation trajectory method [14] by in- 
troducing the fluctuations of the trajectories in the 
subensembles. 



II. THE PROBLEM 

Particle motion in a 2-dimensional stochastic velocity 
field is described by the nonlinear Langevin equation: 



dx{t) 
dt 



v[x(t),f], x(0)=0 



(1) 



where x{t) represents the trajectory in Cartesian coordi- 
nates X = {xi,X2)- The stochastic velocity field v(x, t) is 
divergence-free: V • v(x, t) = and thus its two compo- 
nents vi an V2 can be determined from a stochastic scalar 
field (/)(x, t), as: 



Wi(x, t) 



dxi 



(2) 



where Sij is the antisymmetric tensor {si2 = 1, £21 = 
£11 = £22 =0). In the studies of turbulence of magne- 
tized plasmas, (/)(x, t) is essentially the potential (</> = 
—(j)^/B where (p^{x,t) is the electrostatic potential and 
B is the magnetic field strength) and in fluid turbulence 
<^(x, t)es is the stream function (63 is the unitary vector 
along the axis perpendicular on the plane (.xi, .T2)). 

The potential (/)(x, t) is considered to be a stationary 
and homogeneous Gaussian stochastic field, with zero av- 
erage and given two-point Eulerian correlation function 
(EC) 



(3) 



where (...) denotes the statistical average over the re- 
alizations of 4>{'K,t). The statistical properties of the 
space derivatives of the potential are completely deter- 
mined by those of the potential. They are station- 
ary and homogeneous Gaussian stochastic fields like 
(f){x,t). The two-point Eulerian correlations of the 
derivatives of 0(x, t) are obtained as derivatives of 
the potential EC, E{x,t). We introduce the notation 
-Ei..;fc..(x,i) = {<l)i„{xi,ti) (pk.X'X-i +x,ti +t)) where 



(j>i{-x.,t) = ^0(x, t) and the subscript of E contains the 
indices of the derivatives of the potential in Xi , t\ (left 
factor) separated by semicolon from the indices of the 
derivatives of the potential in xi + x, ti+t (right factor) . 
The absence of indices correspond to a factor (j) inside the 
average (see Eqs. (5) for examples). One obtains 

i^,..;...(x,t) = (-l)"A..._^...^(.,,) (4) 

with n equal to the number of derivatives of the first 

factor (xi, ii) inside the above average. In particular, 
the velocity v(x, t) is such a stationary and homogeneous 
Gaussian stochastic field. The correlation of the velocity 
components and the potential-velocity correlations are 
obtained using the definition (2) of the velocity and Eq. 
(4) as 



(t;i(O,0)t;i(x,t)) = E2;2{^,t) = -^E{x,t), (5) 
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{V2{0,0)V2{x,t)) = £l;l(x,i) = —^^E{x,t), 



(wi(0,0)W2(x,t)) =-^2;l(x,t) = 
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dx\dx2 
d 



^^^■^ ^''dx. 



E{x,t), 
-Eix,t), 



d 



{(l){O,0)vi{x,t)) =SijE.j{x,t) = £ij—E{x,t). 

These correlations (4)- (5) will be used in the following 
calculations. 

The potential is a continuous function of x and t in 
each realization and it determines an unique trajectory 
as the solution Eq. (1). Starting from the above statis- 
tical description of the stochastic potential and from an 
explicit EC, E{x, t), one has to determine the statistical 
properties of the trajectories. The later can be obtained 
from the Lagrangian velocity correlation (LVC), defined 
by: 



Lij{t) = {v,[x{0),0])v,[x{t),t]) 



(6) 



for a stationary process. The mean square displacement 
(a;?(i)) and its derivative, the running diffusion coeffi- 
cient Di{t), are determined by this function [21] as: 



<a;2(t)> = 2 fdrLuir) (t-r), 
Jo 

Di{t) = I dr Lu{t). 
Jo 



(7) 



(8) 



The probability distribution function of the trajectories 
can be obtained from the equation for particle density, 
once the LVC is known [22]. Thus, the solution of the 
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above Langcvin problem consists essentially in determin- 
ing the LVC (6) corresponding to a given EC (3) of the 
stochastic potential. 

This kind of Langevin problem, sometimes named in 
the literature diffusion by continuous movements, is non- 
linear due to the space dependence of the potential, which 
leads to x-dependence of the EC (3). The importance of 
the nonlinearity is characterized by the Kubo number 
defined by 



where V is the amplitude of the stochastic velocity, Tc 
is the correlation time and Ac is the correlation length. 
These parameters appear in the EC of the velocity as 
the maximum value in the origin [V'^ = £^^^(0,0)] and 
the characteristic decay time and length of this functions. 
The Kubo number is thus the ratio of Tc to the average 
time of flight of the particles over the correlation length, 
Tfi = Xc/V. It measures the particle's capacity of ex- 
ploring the space structure of the stochastic velocity field 
before it changes. 

For small Kubo numbers the time variation of the ve- 
locity field is fast and the particles cannot "see" the space 
structure of the velocity field. The condition K <^ 1 
(to -C Tfi) defines the quasilinear regime (or the weak 
turbulence case) for which the results are well estab- 
lished: the diffusion coefficient is Dqi = ( }^/tc)K'^ and 
the trajectories have Gaussian distribution. 

For K > 1 {tc> Tfi) the time variation of the stochas- 
tic potential is slow and the trajectories approximately 
follow the contour lines of 0(x, f). This produces a trap- 
ping eff'ect : the trajectories are confined for long periods 
in small regions. A typical trajectory shows an alterna- 
tion of large displacements and trapping events. The lat- 
ter appear when the particles are close to the maxima or 
minima of the potential and consist of trajectory winding 
on almost closed small size paths. The large displace- 
ments are produced when the trajectories are at small 
absolute values of the potential. Thus there is a strong 
infiuence of trapping on individual trajectories. Trajec- 
tory trapping appears for ii" > 1 and becomes stronger 
as K increases up to the limit of static field {K, Tc = oo) 
where the trapping is permanent. It determines the de- 
crease of the diffusion coefficient and the change of its 
dependence on the Kubo number from the Bohm scal- 
ing [23], [24], Db ^ { Xc/tc)K, to a trapping scaling, 
Dtr ^ ( XI/tc)K^ with 7 < 1. This process was sta- 
tistically described in [14], [15] and the exponent 7 was 
evaluated for given EC of the potential. In the limit of 
static potential field {K, Tc = 00) the transport is subd- 
iffusive with D(t) and {x'^(t)) ~ with a < 1. It 
can be shown (see Section V) that a = 7 and thus the 
same power law describes the large time behavior of the 
mean square displacement in a static potential and the 



dependence on K of the asymptotic diffusion coefficient 
in a time dependent potential. 

Besides this influence on individual trajectories, the 
trapping has collective effects and it generates trajectory 
structures similar with fluid vortices. The trapping ap- 
pears coherently for bundles of neighboring trajectories 
leading to eddying regions. We analyze here the statis- 
tical characteristics of these trajectory structures. We 
show that the dispersion of the trapped trajectories sat- 
urates and that the mean square of the distance between 
trajectories evolve slowly and eventually saturates. The 
motion of trapped particles is almost coherent and leads 
to structures which do not contribute to the transport 
(trajectory dispersion) nor to the relative transport. 

III. THE NESTED SUBENSEMBLE METHOD 

We start from the main idea of the decorrelation tra- 
jectory method [14] . It consists in studying the Langevin 
equation (1) in subensembles (S) of realizations of the 
stochastic held, which are determined by given values of 
the potential and of the velocity in the starting point of 
the trajectories: 

(S): .^(0,0) = <^o, v(0,0)=vO. (10) 

We note that similar subensemble averages of the Eu- 
lerian stochastic velocity field were studied in Ref. [25] 
with the aim of showing that eddies and structures exist 
even in isotropic turbulence. Subensemble Lagrangian 
averages are estimated in [26] on the basis of a rather 
strong assumption and in [12] for a model of rotating 
fluid layers. Our approach is different. 

The statistical properties of the stochastic potential 
and velocity, reduced in the subensemble (S) defined by 
condition (10) are derived in [14]. They are Gaussian 
fields but non-stationary and non- homogeneous, with 
space and time dependent averages and dispersions. The 
averages depend on the parameters of the subensemble 
and are defined by 

cl.^(x,i;5) = (0(x,t))s, V^{x,t;S) = {v[x,t])s (11) 

where denotes the average taken on the realiza- 

tions in (S) and the superscript E is used to underline 
the Eulerian nature of these quantities. They are equal to 
the corresponding imposed condition (10) in x = and 
t = and decay to zero at large distance and/or time. 
The mean squares of the potential and velocity fiuctu- 
ations are zero in x = 0, t = and increase up to the 
values corresponding to the whole set of realizations at 
large distance and/or time. The existence of an average 
Eulerian velocity in the subensemble determines an aver- 
age motion (decorrelation trajectory). This decorrelation 
trajectory is estimated in each subensemble and an ap- 
proximation for the LVC is derived in terms of the set of 
these trajectories. 
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More precisely, the LVC (6) for the whole set of re- 
alizations is obtained by summing up the contributions 
of each subensemble. The latter can be written as 



(«,[x(0),0] Vj[Mt),tl 



vf {vj [x(t),i])„. Thus, the 



problem of evaluating the LVC reduces to the determina- 
tion of the average Lagrangian velocity in each subensem- 
ble (S), V^{t-S) = {-v[yi{f)J])g. This is one of the ad- 
vantages brought by the subensemble analysis: the two- 
point LVC (6) can be expressed as a function of one-point 
averages (V-^(t; 5*)) in subensembles which correspond to 
given initial velocity: 

Lij {t) = j j dct>'' rfvO Pi (S) v° (i; S) (12) 

where Pi{S) is the one-point Gaussian probabil- 
ity density for the values of the potential and ve- 
locity in the point (0,0), defined by Pi{S) = 
{6 [(/)° - (/)(0, 0)] S [vO - v(0, 0)] ) . It represents the 
probability that a realization belongs to the subensemble 
(S) and is obtained as: 

PiiS) 



X exp I — 



2E{0,0) 2£;i;i(0,0) 2£;2;2(0,0) 



(13) 



because the potential and the velocity components are 
not correlated in the same point (((/)(0, 0)wi(0, 0) ) = 
when -E(x) has a maximum in x = 0, as happens in most 
cases). Eq. (12) is an exact equation. 

The approximation introduced in the decorrelation tra- 
jectory method is in the estimation of 'V^{t;S) and es- 
sentially consists in neglecting the fluctuations of the tra- 
jectories around the average trajectory in (S). Thus the 
average Lagrangian velocity in (S) is approximated with 
the average Eulerian velocity calculated along the aver- 
age trajectory: 



S ' 



I.e. 



V^{t;S)^V^ [X(i,S'),t; 5*] 



(14) 



(15) 



where X(f, S) is the average trajectory in (S), X(t, S) = 
{x{t))g . Then, this average trajectory in (S) (decorrela- 
tion trajectory) is determined from the equation 



rfX(t, S) 
Jt 



= V^ [X(t,5),i;5] 



so that an explicit expression for the average Lagrangian 
velocity is obtained. It was shown that Eq. (16) is of 
Hamiltonian type with the Hamiltonian function equal to 
the average Lagrangian potential (f>^ [X(t, S),t; S] . Thus, 
in the static case the average trajectories obtained from 



the approximation (14) are periodic and evolve on closed 
paths. They provide a representation of the trapping 
and lead to a subdiffusive transport for the static poten- 
tial and to trapping scaling for the asymptotic diffusion 
coefficient, Dt^ ~ K'' with 7 < 1, in time dependent 
stochastic potentials with slow variation {K > 1). 

The approximation (14) seems to be rather rough but, 
because it is performed in the subensemble, there are 
two aspects which contribute to improving its accuracy. 
One is due to the fact that the fluctuations of the ve- 
locity in (S), which are the source of the trajectory fluc- 
tuations, are zero in the starting point of the trajecto- 
ries and become important only if the trajectory reaches 
large enough distances. The second is related to the fact 
that the trajectories in the subensemble are superdeter- 
mined. Besides the necessary and sufficient initial condi- 
tion x(0) = 0, they have supplementary initial conditions 
determined by the definition (10) of the subensembles. 
This reduces the differences between the trajectories in 
(S) and thus the fluctuations. The first description of 
the trapping process in qualitative agreement with the 
numerical simulations was obtained using this approxi- 
mations [14]- [18]. However, there are important qualita- 
tive aspects that are not obtained from this approxima- 
tion. The most evident concerns the average trajectory 
in the subensemble. In the static case the trajectory in 
each realization (solution of (1)) is periodic but the aver- 
age of such trajectories cannot be periodic (as obtained 
from (16)) since they have different periods (distributed 
around some average value which depends on (S)). An- 
other aspect is discussed in [19], [20] and concerns the 
average Lagrangian velocity in a biased stochastic poten- 
tial. It is thus necessary to improve the decorrelation tra- 
jectory method by taking into account the fluctuations of 
the trajectories in the subensembles. This development 
is also required by the aim of this paper. The decorrela- 
tion trajectory method as presented in [14] is not able to 
describe the statistics of the trajectories in (S) nor the 
two-particle statistical quantities. 

The analysis of the decorrelation trajectory method 
leads to the conclusion that this method succeeds in de- 
scribing the trajectory trapping due to the fact that it 
maintains the invariance of the average Lagrangian po- 
tential in (S). The motion described by Eqs. (1), (2) has 
the velocity at any moment perpendicular to the gra- 
dient of the potential in the instantaneous position x(i) 
and the time variation of the Lagrangian potential is pro- 
duced only by the explicit time dependence of 0(x, t) 



(16) d<i)[^{t),t] 



dt 



d(j)[yi{t)] dxi _^ d(l)[x{t),t] 



d<P[x{t),t] 



dxi dt 



dt 



dt 



(17) 



Thus, for the static case where d<p/dt = (tc, K 00), 
the potential is an invariant of the motion. The trajec- 
tories are on the contour lines of the potential and the 
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motion is periodic. For slowly varying or large ampli- 
tude potentials such that Tc > r/; (i.e. K > 1), dcp/dt 
is small and there is an approximate invariance of the 
potential along the trajectory in each realization. Tra- 
jectory trapping is essentially related to this invariance 
of the Lagrangian potential. 

The average Lagrangian potential in (S) is invariant in 
the frame of the decorrelation trajectory method. Indeed, 
one obtains, by neglecting trajectory fluctuations as in 
(14), the average Lagrangian potential in (S) as 

$^(i; S) = {<j> [x(t), t]) ^ [X(t, S), t; S] . (18) 

Since the Eulerian quantities V^(x, t; S) and ^^{x,t; S) 
are related by an equation similar to (2) 



y,^(x,i;S)=£, 



^^^(x,^;^) 

' dxi ' 



(19) 



the time derivative of (18) is zero for the static case 
and, due to the definition (10) of the subensembles (S), 
$^(t; S) = $^(0; S) = (fP at any time. 

The aim of this paper is to consider the fluctuations 
of the trajectories in the subensemble (S) and to study 
their effect. The approximations (14), (14) that neglect 
trajectory fluctuations have to be replaced. As suggested 
by the above discussions, this development must be per- 
formed having in mind the idea of maintaining the con- 
dition of invariance of $^(t; S). Moreover, the invariance 
of the Lagrangian potential applies to each trajectory, 
i.e. in each realization. Thus, it determines, besides 
the invariance of $^(t;5), other statistical constraints. 
Actually, in the static case the probability distribution 
function for the Lagrangian potential in a subensemble 
(S) is 



,<)= (n^--^ [x(t)]})s = (5(0-0°) 



(20) 



because [x(t)] = for all trajectories in (S) and for any 
time moment. We note that the knowledge of P^{4>, t) is 
not very helpful in solving the problem of determining the 
statistical properties of the trajectories but it rather im- 
poses strong conditions on the possible approximations. 
For example, the hypothesis that trajectory fluctuations 
in (S) are Gaussian around the average trajectory is ex- 
cluded by Eq. (20). It can be shown that such distri- 
bution is not compatible with the ^-distribution of the 
Lagrangian potential. 

Trajectory fluctuations in the subensembles (S) are 
considered here in agreement with the condition (20) by 
separating the realizations in (S) in subensembles (S2) 
corresponding to fixed values of the second derivatives of 
the potential in x = 0, t = 



(52): 0i,(O,O) 



dxidxj 



(21) 



x=O,t=0 



where ij = 11, 12, 22. The Langevin equation (1) is stud- 
ied in these subensembles (S2). The average trajectory is 



approximately determined by neglecting trajectory fluc- 
tuations in (S2). Pushed to the subensembles (S2), this 
approximation is much more accurate than taken in (S) 
because the trajectories in (S2) are much more superde- 
termined than those in (S): three supplementary initial 
conditions (21) are added to the initial conditions (10). 
Moreover, the source of trajectory fluctuations, the ve- 
locity fluctuations, has smaller amplitude in (S2) than in 
(S). Thus the accuracy of this method is much increased 
compared to the accuracy of the decorrelation trajectory 
method. But the main advantage of performing this de- 
velopment consists in the fact that it provides for each 
subensemble (S) a collection of trajectories instead of one 
decorrelation trajectory. It is thus possible to determine 
the statistical properties of the trajectories in (S) by per- 
forming averages over the subensembles (S2) contained 
in (S). We determine using this nested subensemble pro- 
cedure the dispersion of the trajectories in (S) and their 
probability distribution function, besides the average tra- 
jectory. The statistical evolution of the distance between 
two trajectories is also studied by this method and the 
average, the dispersion and the probability distribution 
function are determined as functions of time. 

We note that this nested subensemble approach can be 
further developed by introducing subensembles deflned 
by higher order derivatives of the potential, (S3), (S4), 
.... This systematic expansion fulflls at each order higher 
than 1 all the conditions required by the invariance of 
the Lagrangian potential. It is however expected that the 
main statistical properties of the stochastic trajectories in 
(S) are already obtained at the second order and that the 
higher orders contribute with corrections to these results. 
The nested subensembles (S), (S2) are considered in this 
study. 

The explicit calculations consist of the following steps. 
First, the statistical properties of the stochastic potential 
and velocity, reduced in the subensemble (S2) defined by 
conditions (21) and (10) are derived. Namely, the prob- 
ability that a realization belongs to a subensemble (S2) 
and the subensemble average Eulerian velocity and po- 
tential are determined. These average quantities lead to 
the equation for the average trajectory in (S2) by neglect- 
ing the fiuctuations of the trajectories in (S2). Then the 
statistical properties of the trajectories and of the dis- 
tance between two trajectories in the subensemble (S) 
are expressed as functions of the average trajectories in 
all subensembles (S2) contained in (S). The calculations 
are for the static stochastic potential corresponding to 
the strongest trapping. The transport in time-dependent 
potential is discussed in Section VI B. 



A. Eulerian statistics in the subensemble (S2) 

The (one-point) probability that a realization of the 
stochastic potential belongs to the subensemble (S2) con- 
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tained in the subensemble (S) is defined by 



Pi{S2) = 



(22) 



where the product is for ij = 11, 12, 22. It is calculated 
using the Fourier representation of the 5-functions and 
performing the average of the resulting exponential of the 
stochastic Gaussian quantities (see e.g. [27]). One obtains 
after straightforward calculations: 



Pi{S2) = 



1 



(27r) 



3/2 



12:12 



exp 



[E 



2E, 



-1/2 



2c 



(23) 



gjci aia2Ci2 
2c c 



where Cj are constants given by 

Ci = £n;n(0)-SL(0)/£;(0), 
C2 = E22-MO)-El2;iO)/E{0), 

C12 = Sll;22(0) - Ell.{0)E22;{0)/E{0) 



(24) 



C = C1C2 



c 



12 



and tti are essentially the parameters of the subensem- 
ble (S2) 



ai 
a2 



^?i 
4 



iPE^i,{Q)/E{0), 
fE22MIE{0). 



(25) 



The average Eulerian potential in the subensemble 
(S2), <I>^(x; 52) = {(f>{x))g2 , is determined by the condi- 
tional average corresponding to (21): 



$^(x; S2) 



(26) 



(0(x) S - ^(0)] ^ [yo - v(0)] n <5 [^°, - <l>ij (0)] ) 
Pi{S)Pi{S2) 

This average is calculated by introducing the Fourier rep- 
resentation of the ^— functions which leads to the average 
of (/'(x) multiplied with an exponential of a linear combi- 
nation of 4>{0), v(0) and 0ij(O). This average is obtained 
as the derivative at a parameter a taken in a = of the 
average of the exponential of the above linear combina- 
tion with an additional term a(/)(x). One obtains after 
performing the inverse Fourier transforms 



$^(x; 52) 



Ei^)^ + Ei,{^)^+Eij.{^) 



(27) 



Pi(5)Pi(52) 



Pi(5)Pi(52) 
which can be written explicitly as 



(28) 



ai {E22;{0)C12 - Eu;{0)C2) _^ 02 (^^11; (0)ci2 - ^^22; (0)ci ) 



v'2 



^12; (X) ^0^^ 



i^2;2(0) ' Ev,l{0) '^^12;12(0) 



■£;ii;(x) (aiC2 - a2Ci2) _^ ^22; (x) (a2Ci - aici2) 
c c 

The subensemble (S2) average potential (28) equals 0° 
in X = 0, t = and it decays to zero at large x, as the 
average Eulerian potential in the upper subensemble (S). 
Its expression is more complicated and depends on the 
second order derivatives cf)'^j that label (S2). The poten- 
tial considered only in the realizations in the subensemble 
(S2) is a non-stationary and non-homogeneous Gaussian 
field having a space-time dependent average. 

As in the whole ensemble and as in (S), the statistical 
properties of the velocity field in the subensemble (S2) 
are deduced from those of the potential in (S2). The 
velocity in the subensemble (S2) (21) is a non-stationary 
and non-homogeneous Gaussian stochastic field having 
a space-time dependent average. This average Eulerian 
velocity is calculated by the same procedure used for the 
potential (28). The relation (2) between velocity and 
potential in each realization holds between the respective 
average quantities calculated in the nested subensembles. 
It was obtained in the subensemble (S), Eq. (19), and it 
can be shown that : 



y,^(x; 52) = e, 



9$^(x; 52) 



(29) 



Thus the average velocity in the subensemble (S2) is 
divergence-free: V • Vf (x; 52) = 0. 

It is interesting to note that the potential and the ve- 
locity in the subensembles (S) and (S2) arc deterministic 



quantities in x = (0(0) 



v(0) 



v" for all re- 



alizations in (S), thus also in (S2)). As |x| grows, the 
average values decay to zero and the fluctuations build 
up progressively and eventually become the same as in 
the global statistical ensemble. The amplitude of the 
fluctuations in a point x is smaller in the subensemble 
(S2) than in the subensemble (S). 

This nested subensemble procedure evidences, in the 
zero-average stochastic velocity field, a set of average 
velocities (corresponding to each subensemble). They 
depend on the statistical characteristics of the veloc- 
ity field (the correlation and the constraint imposed in 
the problem, i.e. the zero- divergence condition). In the 
nested subensembles introduced here, the following re- 
lation holds between the (S2) average velocities and the 
(S) average velocity: 
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V^(x;5)= j dcj>\^dct>l^d<t>l2Pi{S2)\^{^;S2). (30) 

Similar equations can be written for all statistical quan- 
tities defined in the nested subensembles. 



B. Average Lagrangian velocity in the subensemble 

(S2) 

The average Eulerian velocity (29) determines an av- 
erage displacement in the subensemble (S2), X(t; 5*2). It 
is the time integral of the average Lagrangian velocity in 
(S2), V^(t; 5*2). The latter is evaluated using an approx- 
imation similar to (14) as 



V^(i; 52) ^ V^[X(i; 52); 52] 



(31) 



by neglecting the fluctuations of the trajectories in the 
subensemble (S2). A nonlinear equation for X(t; 52) is 
so obtained 



o;X(^, 52) ^ [X(t;52);52]. 



(32) 



With this approximation the average Lagrangian poten- 
tial in (S2) is {if [yi{t)])s2 - [Mt; 82); 52] and due to 
Eq.(29) 

[X(t; 52); 52] = ^u^*'' [X(t; 52); 52] . 

which shows that Eq. (32) has a Hamiltonian struc- 
ture. The solution of this equation with the initial 
condition X(0; 52) = ensures the invariance of its 
time-independent Hamiltonian function <I>^ [X(i; 52); 5]. 
Since the Eulerian average potential (28) has the value 
(/>° in X = 0, the average Lagrangian potential equals (fP 
at any time moment. 



C. Lagrangian statistics in the subensemble (S) 

The nested subensemble method provides for a 
subensemble (S) an ensemble of trajectories, the aver- 
age trajectories X(t; 52), one for each subensemble (S2) 
contained in (S). The statistical properties of these tra- 
jectories can be determined performing averages over the 
subensembles (S2). 

The average velocity in (S), 'V^{t;S), is determined 
by averaging V^(t;52) obtained from the solution of 
Eq.(32) over all subensembles (S2) contained in (S) 

V^(i; 5) = y #?i#;2#22A(52)V^ [X(t; 52); 52] . 

(33) 



Similar equations hold for all the statistical quantities 
corresponding to the subensemble (S). The average tra- 
jectory in the subensemble (S), X(f; 5) = (x(f))g , is de- 
termined by averaging X(f; 52), the solution of Eq. (32) 



X(t;5)= / d(/)°i^d(f>"^2d<p22PiiS'2)Mt;S2). 



(34) 



The dispersion of the trajectories, 
di{t;S) = (^{xi{t)-X{t;S)f)^ is 

di{t; S) = j #;i#?2#02A(52)X2(t; 52) - X^{t; 5). 

(35) 

The probability distribution function (pdf) for the tra- 
jectories in (S) is determined by integrating the pdf in 
the subensemble (S2) which in this approximation is 
(5 [x - X(t; 52)] 

P^(x,t) = y#;i#?2#02-Pi(52)5[x-X(t;52)]. (36) 

The pdf for the Lagrangian potential <p [^(i)] i^i the 
subensemble (S) obtained by this method equals 5{4>—(jP) 
since the potential is equal to (fP on all trajectories con- 
sidered in the average. This shows that the approxima- 
tion (31) introduced in this method ensures entirely the 
statistical properties of the Lagrangian potential. 

The statistical properties in the subensemble (S) of 
the distance between two trajectories 5{t) = x'(t) — x{t) 
starting from x'(0) = 5q and x(0)= respectively can 
also be determined using the average over the subensem- 
bles (S2). The average trajectory in (S2), (x'(t))g2, 
is the solution of Eq. (32) with the initial condition 
X'(0;52) = 6o. The average, the dispersion and the pdf 
of S{t) in the subensemble (S) are determined by equa- 
tions similar to (34)-(36) where X(t; 52) is replaced by 
X.'{t;S2)-X{t;S2). 



D. Running diffusion coefficient 

The correlation of the Lagrangian velocity in the whole 
set of realizations is determined using Eq. (12) where 
V^(t; 5) is the time derivative of the average trajectory 
in (S), given by Eq. (34). The running diffusion coeffi- 
cient (8) is obtained from (12) and (34) as 

Di {t) = j j dcjP dv° Pi (5) v° Xi{t;S). (37) 

In the case of an isotropic stochastic field, the integral 

over the orientation of the velocity v" can be performed 
analytically [14] and one obtains for the static case 



m = 



1 



2tt ^£(O,0)£;i;i(O,0) 



(38) 
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and 

m = D'{t) , (39) 

where Xi{t; S) is the component of the average trajec- 
tory along v°, determined from Eq.(34), and D'{t) is the 
derivative of the function D{t). We note that the same 
analytical expression for D{t) in terms of the average 
trajectory Xi{t;S) is obtained in [14] by means of the 
decorrelation trajectory method. But the average trajec- 
tory was determined there as solution of a Hamiltonian 
equation while here it is the average (34) of the average 
trajectories in (S2). 

IV. EXPLICIT CALCULATIONS 

The nested subensemble method actually is based on 
the classification of the stochastic trajectories in groups 
(subensembles) according to some resemblance between 
them. The most important criterion in this classifica- 
tion is the value of the potential in the starting point 
of the trajectories. All trajectories contained in such a 
group evolve on contour lines with the same value of the 
potential. Consequently their paths and periods are sta- 
tistically similar in the sense that they have an average 
size and period. Thus this condition determines a global 
resemblance of the trajectories (extended at long time). 
This condition is imposed beginning with the first level 
of classification (in the subensembles (S)). Other criteria 
of the classification are the velocity and the derivatives of 
the velocity in the origin. These are not conserved quan- 
tities ant they influence the shape of the trajectory only 
at small time for time intervals that grow with the num- 
ber of imposed conditions. The value of the initial veloc- 
ity is fixed in the subensemble (S), then each subensem- 
ble (S) is divided in smaller subensembles (S2) according 
to the value of the derivatives of th(^ v(^locity (second 
derivatives of the potential) . This classification can con- 
tinue in principle and at each step the resemblance of 
the trajectories contained in a group is increased and the 
number of groups grows. The approximation consists in 
neglecting the differences between the trajectories in a 
group. With this condition it is possible to determine 
an average trajectory for each subensemble. Thus the 
nested subensemble method determines a set of trajecto- 
ries X(t; S2) and a weighting factor for each one. Then, 
the statistical properties of the stochastic trajectory are 
obtained by performing averages over these trajectories. 

Except for some special case, the trajectories X(t; 5*2) 
have to be numerically calculated by solving the Hamil- 
tonian system (32). This procedure appears to be very 
similar with a direct numerical study of the simulated 
trajectories. There are however essential differences. The 



average trajectories are obtained from a rather smooth 
and simple time-independent Hamiltonian. They are pe- 
riodic functions and thus are calculated only for a pe- 
riod. The number of trajectories is much smaller than 
in the numerical study due to the weighting factor deter- 
mined analytically. This reduced very much the calcula- 
tion time, such that it can be performed on PC. Moreover 
such a calculation performed for a static stochastic po- 
tential with a given EC determines the solution for the 
time dependent potential with arbitrary time factor in 
the EC (see Section VI B). 

We have developed an algorithm for calculating the 
statistical characteristics of the trajectories in subensem- 
bles (S) and the running diffusion coefficient (38) for 
given EC of the potential. The trajectories are calculated 
for a period using a variable integration step determined 
monitored by the precision obtained for the potential. 
Values of this precision of 10~^ — 10~^ ensure the stability 
of the calculated D{t). The order of performing the in- 
tegrals in Eq.(38) appears to be important. The integral 
over u is first calculated. This parameter is factorized in 
the expression of the average Eulerian potential in (S2) 
(28) such that it appears as a factor in the right hand 
side of Eq.(32). This factor is introduced in the time 
variable and so only the trajectories with u = 1 need 
to be calculated. The values of the function Xi{ut;S2) 
are determined by interpolation and using the periodic- 
ity, the range of u is determined from the range of the 
exponential factor and the step du is determined at each 
integration such that a large enough number of points 
(30-50) exists on each period. When there are more than 
about 50 periods on the range of u the integration is not 
performed because its value is negligible. The next in- 
tegrations are over (I)i2, and (pl2- The integrand for 
each of these integrals is over the range determined by 
the exponential factor and we have calculated it on a 
mesh with constant step with 31-61 points. The last is 
performed the integral over (jp. Due to trajectory trap- 
ping, the range of this integral is continuously reduced 
as time increases (ct>max ^ when t — > oo). The range 
of this integral is calculated clS 3j function of time for 
the interval of calculation of D{t). The calculations start 
with a large value of '/'maxi obtained from the expo- 
nential factor. At the time when the integrand becomes 
approximately zero on half of this range f/'^ax is reduced 
and the integration of the trajectories is taken again from 
t ~ for the new values of (fP . The mesh for (fp' has vari- 
able steps that increase toward large values of because 
the function has strong variations at small 0°. The tests 
performed with this code have shown that the numerical 
calculations are rather fast and accurate and they can be 
advanced up to large values of time. For instance, using 
the decorrelation trajectory method the duration of the 
calculation of D{t) up to time of the order 10~2 is of the 
order of 10 seconds on an usual PC. Using the nested 
subensemble method the calculation time is of the order 
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of one day because the number of calculated trajectories 
increases with a factor 10"*. 

The nested subensemble method determines the LVC 
for test particles moving in a stochastic potential with 
given EC. The main condition for using this method 
is that the transport is stationary, which usually cor- 
responds to stationary and homogeneous stochastic po- 
tentials. The potential field is Gaussian. The time- 
dependent diffusion coefficient (38) corresponds to an 
isotropic potential but this is not a restriction for this 
method. The above calculation are for a static potential 
but they can be extended to time dependent case (see 
Section VI B and [14]). 



V. TRAJECTORY STRUCTURES 

We present in this Section typical results obtained 

for the statistical characteristics of the trajectories in a 
subensemble (S2). We need to specify the EC of the 
stochastic potential, which we choose as 



£;(x) 



l + {xl + xl)/2- 



(40) 



This is the EC of a normalized stochastic potential with 
amplitude E{Q) = 1, Ac = 1 and the time with Tfi — 1. 
The velocity v° that defines the subcnscmbles (S) is taken 
along the x\ axis. 

The average trajectory in the subensemblcs (S2), solu- 
tion of Eq.(32), is a periodic function of time and evolves 
on a closed path for most of the subcnscmbles (S2). 
There arc also some open paths for subcnscmbles with 
(f)"^ — Q and trajectories with very large periods for small 
values of \<t>^\- Some examples of paths of the average 
trajectories in (S2) arc presented in Fig. 1. There is a 
clear difference between the trajectories corresponding to 
small |(/)°| (Fig. la for = 0) and large \(jP\ (Fig. lb for 
(ffi = 1). In the first case there are open trajectories, large 
displacements and large periods for the closed paths. In 
the second case the trajectories have small size and their 
periods are much smaller. The size of the path and the 
period of the trajectory depend on the six parameters 
that define the nested subensembles (S), (S2). 

The average trajectory in the upper subensemble (S) 
is obtained from (34). Typical average trajectories in 
subensemble (S) are presented in Fig. 2. The time depen- 
dence of X(t; S) is presented in Fig. 3. These trajectories 
are not periodic. They evolve on spiral shaped paths, ex- 
cept for the subensemble with 0° = which yields a 
continuous displacement along v°. The size of the paths 
depends on the parameters of (S), (jP and u = |v°| : it is 
large for small | (jP | and large u and it decreases as | 
increases. The displacement along the initial velocity v 
decays to zero in a characteristic time while the dis- 
placement perpendicular to v° saturates at a finite value 



whose sign is the same as the sign of (fp . The saturation 
time Ts depends on the parameters of the subensemble 
(S). It increases when \(jp\ decreases (as the size of the 
paths) and when (fP ^ Q \t becomes infinite. Thus the 
average trajectories in (S) obtained here are completely 
diff'erent from those obtained by means of the decorrela- 
tion trajectory method [14]. The later arc periodic fimc- 
tions of time and their paths are closed (see Fig. 2 for 
comparison). This means that the fluctuations of the tra- 
jectories have a strong influence on the average trajectory 
in (S). They determine the time-saturation of the aver- 
age trajectory in (S) by the mixing of the closed periodic 
trajectories. 

The dispersion of the trajectories in the subensemble 
(S) obtained from Eq. (35) is presented in Fig. 3 as a 
function of time for 0° = (Fig. 3a) and (fp = \ (Fig. 
3b). One can see that in the first case the dispersion con- 
tinuously increases (Fig. 3a) while in the second case it 
saturates after a more complicated evolution (Fig. 3b). 
The saturation time is the same as for the average tra- 
jectory. The amplitude of the trajectory fluctuations is 
comparable with the average displacement. 

The trajectories in the subensemble (S) are Gaussian 
at small time t <^ Tfi but their distribution is strongly 
modified as time increases. The pdf obtained from Eq. 
(36) is represented in Fig. 4. Important differences can 
be observed between the subensembles with (jP = (Fig. 
4a) and those with large 10*^1 (Fig. 4b). In the first case 
the pdf is symmetric around v'', it develops a narrow 
maximum in x = and an annulus that expands contin- 
uously in the direction of v° as time increases (Fig. 4a). 
The velocity of this part of trajectories is larger than the 
average velocity. The path of the average displacement is 
also represented in Fig. 4a. The end point of this curve 
is the average position at the moment corresponding to 
the representation of the pdf (t — lOOr/;). It is located 
between the two maxima of the pdf in a region where the 
later is practically zero. The pdf for subensembles with 
large |0°| is completely different. It saturates in a time Ts 
at a function that has a narrow maximiim in x = and 
extends only on a small region (with 2:2 > for (fp > 0). 

Thus the statistical characteristics of the trajectories 
in a subensemble (S) with a large values of | | are com- 
pletely different of those corresponding to subensembles 
with (jP = 0. The average, the dispersion and the pdf 
of these trajectories saturate. This shows that there is 
a quasi-coherent motion in such subensembles consist- 
ing in trajectory rotation on closed paths, with localized 
pdf and small saturated dispersion. The trajectories form 
structures similar with fiuid vortices. These structures or 
eddying regions are permanent in static stochastic poten- 
tials. The saturation time Tg represents the average time 
necessary for the formation of the structure. The aver- 
age size of the structure is represented by the asymptotic 
average displacement |X(t;S')| a,t t ^ Tg. The disper- 
sion of the trajectories in the structure is given by the 
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asymptotic value of di{t;S). We have found that these 
characteristic parameters of the trajectory structures de- 
pend on the parameters of the subensemble (S). The size, 
the dispersion and the build up time of the structures in- 
crease when |(^°| decreases and go to infinity at (/)° = 0. 

The existence of the trajectory structures is confirmed 
by the statistical properties of the distance between two 
neighboring trajectories, 6{t) = x'(i) — x(i). Typical re- 
sults obtained for the second moment of the relative dis- 
placement 5{t), arc presented in Fig. 5 (con- 
tinuous lines) compared with the second moments of the 
absolute displacements, (dashed lines). In the 
subcnscmbles with large l^*^! (Fig. 5b). the evolution 
of (^S^{t))g shows that it maintains long time the ini- 
tial value and that it reaches values comparable with 
the absolute displacement (a;|(t))g only after a very long 
time (of the order of lOOr/;). Thus the relative motion 
is strongly hindered and a very strong clump effect ap- 
pears in the subensembles with large • There is a 
very strong degree of coherence of the relative motion for 
these trajectories showing that they form structures. On 
the contrary, the clump effect is very weak (practically 
absent) for the trajectories which are not in the struc- 
tures (those in the subensembles with = 0). As seen 
in Fig. 5a, (6f{t))g has values comparable to those of 
(a;f(t))^ since the first stage of the evolution, at time 
much smaller r//. The pdf of S(t) in a subensemble (S) 
shows that the relative motion is not Gaussian. In the 
case of structures (large the pdf remains very lo- 
calized around zero and saturates (Fig. 6b). In the case 
of free trajectories {(p^ = 0), the pdf has a more compli- 
cated shape and extends continuously. At large times it 
is similar with the pdf of the trajectories (Fig. 6a). 

Thus, the trapping of the trajectories has a strong in- 
fluence on the statistical characteristics of the relative 
motion. It produces an anomalous clump effect. In 
the absence of trapping the clump effect appears as an 
exponential time dependence of the average square dis- 
tance between two trajectories, of the type {Sf{t)) = 
Sq exp{t / Tci) J where the clump characteristic time Td is 
of the order of the diffusion time which is the flight time 
Tfi at large K (see [28] and the references therein or the 
review paper [4]). The distance between two neighbor- 
ing trajectories remains close to the initial value during 
a time Td and then ((5f(t)) grows rapidly and reaches 
a diffusive behavior with the diffusion coefficient equal 
to 2D. The process of trajectory trapping determines a 
complete change of the clump effect. It appears only for 
a part of the trajectories, those contained in subcnscm- 
bles (S) with large \(l>°\ , and is very strong. The life time 
of the clump is much larger than t/j and than the satu- 
ration time Tg. The time evolution of the relative square 
displacement is much slower. Neighboring particles have 
thus a coherent motion for a long time. For the other 
part of the trajectories, those that move along contour 



lines of the potential with cj)'^ = 0, the clump effect is ab- 
sent and the relative motion become rapidly incoherent, 
after a time interval smaller than r/;. 

VI. TRANSPORT 

A. Static stochastic potential 

The Lagrangian velocity correlation and the time de- 
pendent diffusion coefScient for the whole ensemble of 
trajectories are determined from Eq. (38). The in- 
tegral over the parameters of the subensembles (S) of 
the average displacement in (S) has to be calculated. 
Since Xi{t;S) decays to zero in a time Ts{S), the tra- 
jectory structures have only a transient contribution to 
the running diffusion coefficient. At times larger that 
Ts{S) the contribution of the subensemble (S) vanishes. 
As time increases the diffusion cocfffcicnt D(t) is deter- 
mined by a smaller and smaller number of trajectories, 
those contained in large structures with large saturation 
time. The results obtained for D{t) arc presented in 
Fig. 7 where the dimensionless function F{t) = D{t)/ Db 
is plotted (continuous line). The transport is subdiffu- 
sive in such static stochastic potential. One can observe 
the change that appears at f > r/;. The running diffu- 
sion coefficient begins to decrease and eventually goes 
to zero. A power law decay was obtained at t > Tfi as 
D{t) = V Xc{t/Tfi)~°''^^ . The LVC becomes negative 
at this time and after a minimum it has a long negative 
algebraic tail that decays to zero. The positive and the 
negative parts of L{t) have equal time integral such that 
/o L(r)dT = D(t) 0. The mean square displacement is 
(a;2(t)) ~ and thus the process is subdiffusive. 

The result obtained for D{t) with the decorrelation tra- 
jectory method is also plotted in Fig. 7 (dashed line). It 
is surprisingly close to the result of the nested subensem- 
ble method although the two methods yield completely 
different average trajectories in the subensembles (S) 
(Fig. 2). This shows that by introducing the subensem- 
bles (S2) in the nested subensemble method a strong 
qualitative improvement of the statistical results in the 
next upper subensemble (S) is obtained and only a small 
correction at the level of the whole set of realizations. 
It is thus expected that the development of the method 
by introducing higher order derivatives and the corre- 
sponding nested subensembles (S3), (S4), ...would yield 
only small corrections for the physically interesting re- 
sults that concern the diffusion coefficient D{t) and the 
statistical properties of the trajectory structures. This 
nested subensemble method appears to be fast conver- 
gent. This is a consequence of the fact that the mixing 
of periodic trajectories, which characterizes this nonlin- 
ear stochastic process, is directly described at each order 
of our approach. The results obtained in first order (the 
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decorrelation trajectory method) for D{t) are thus vali- 
dated by the present second order calculations. 

B. Time-dependent stochastic potential 

For a time-dependent potential 0(x, t) (finite Tc and 
K), it is also possible to apply the nested subensemble 
method following the same procedure as above. A very 
simple analytical solution is obtained when the stochastic 
potential has independent time and space variations such 
that its EC is E{x.)h{t). In this case, the average Eulcrian 
potential in the subensemble (S) is given by Eq.(28) mul- 
tiplied with the factor h{t). This factor is transmitted to 
the average Eulerian velocity in (S) (29) and it appears 
in the equation (32) for the average trajectory in (S2). A 
change of variable from t to 

e{t)= I h{T)dT (41) 
Jo 

can be performed in Eq.(32) and thus the average tra- 
jectory in (S2) for a time dependent potential can be 
written in terms of the average trajectory for the static 
case as X(0(i);5'2). The argument 9{t) determined by 
the time-dependence of the potential is 9{t) = t bX small 
t and saturates at a constant which is the decorrelation 
time 6{t) Tc- The same expression (38) is eventually 
obtained for the time dependent diffusion coefficient but 
with Xi{t; S) replaced by Xi{0{t); S) and thus the diffu- 
sion coefficient is 

D'^^it) = D [e{t)] . (42) 

The limit for t ^ c» is finite which shows that the trans- 
port is diffusive in a time dependent stochastic potential 
and the diffusion coefficient is 

D*'^ = D{tc) = DbF{tc). (43) 

This equation shows that the asymptotic diffusion coef- 
ficient is determined by the time dependent diffusion co- 
efficient D{t) corresponding to the static potential {F{t) 
is the function plotted in Fig. 7 and represents the nor- 
malized diffusion coefficient in the static potential and 
Db = {^l/TcjK = VXc is the Bohm diffusion coefficient 
obtained when trajectory trapping is neglected). In the 
limit of small K, the quasilinear result is recovered from 
Eq.(43) and, at large K, I?*'' is reduced compared to the 
Bohm diffusion coefficient by a factor F{K) < 1 which 
accounts for trajectory trapping. For the above EC of the 
potential, Eq.(43) gives the large K scaling law D*'^ « ( 
\1/tc)K^ with 7 = 0.57. The exponent 7 depends on the 
EC of the potential, namely on its space dependence at 
large distances. It is not a fixed value as in the estima- 
tion based on percolation theory [13]. A detailed study 
of the effect of the EC of the potential on the scaling of 
the diffusion coefficient is in progress. 



Thus, the study of the static case permits to deter- 
mine the asymptotic diffusion coefficient in a time de- 
pendent stochastic potential. This property appears in 
the results of the nested subensemble method (and in 
the decorrelation trajectory method) but it is possible to 
demonstrate in general that the time dependence of the 
diffusion coefficient in the subdiffusive static case deter- 
mines the Kubo number dependence of the asymptotic 
diffusion coefficient in a time dependent potential. The 
subdiffusive transport corresponds to Lagrangian corre- 
lations L(t) which have the property: 

D{t) = / L{t)dt 0. (44) 
Jo 

This shows that such a correlation has negative parts 

that compensate the small time correlation which is al- 
ways positive. We suppose that D(t) decays to zero as 
D{t) w {t/Tfi)~" and consequently the LVC behaves as 
L{t) w {t/Tfi)~°'~^ . In the time-dependent case, the vari- 
ation of the stochastic field produces the decorrelation of 
the Lagrangian velocity and consequently the Lagrangian 
correlation decays to zero at t ^ Tc. The asymptotic dif- 
fusion coefficient can be evaluated as 

D= L{t)dt2^ / L{t)dt 
Jo Jo 

and using Eq. (44) one can write 

D ~ - / L{t)dt « (Te/r/O"" = K-'^- 

Jtc 

Thus the exponent a of the time decay of the siibdiffusive 
transport coefficient in the static case determines the ex- 
ponent of the K dependence of the asymptotic diffusion 
coefficient in the time-dependent case. This means that 
Eq.(43) holds even if the evolution of £)*''(<) is not given 
by Eq.(42) as happens for example when the EC of the 
potential is not factorized. 

The time variation of the potential determines a decor- 
relation effect. After a time of the order Tc the config- 
uration of the stochastic potential changes. A competi- 
tion appears between the intrinsic tendency of the tra- 
jectories to form structures and the destruction of these 
structures produced by the time variation of the poten- 
tial field. Structures with Tg {S) > Tc cannot exist and the 
corresponding trajectories produce a diffusive transport. 
Small structures that build up rapidly (with Ts{S) ^ Tc) 
still exist if the correlation time of the field is longer than 
the flight time (tc > Tfi, K > 1). These vortical struc- 
tures do not contribute to the large time values of the 
diffusion coefficient and the transport is reduced. 

VII. CONCLUSIONS 

We have studied the special problem of test parti- 
cle transport in 2-dimensional divergence-free stochas- 
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tic velocity fields, which is characterized by the intrin- 
sic trapping of the trajectories on the contour line of 
the stochastic potential. We have developed a semi- 
analytical statistical approach, the nested subcnscmble 
method. The time dependent diffusion coefficient is de- 
termined by means of a set of deterministic trajectories, 
the average trajectories in subcnscmblcs with given val- 
ues of the potential and of the velocity in their starting 
point. These trajectories are obtained by dividing each 
subcnscmble in a class of subcnscmblcs defined by the 
values of the second derivatives of the potential. Thus, 
the nested subensemble approach reduces the problem of 
determining the statistical behavior of the stochastic tra- 
jectories to the calculation of weighted averages of some 
smooth, deterministic trajectories determined from the 
EC of the stochastic potential. 

The statistical characteristics of subensembles of tra- 
jectories are obtained with this method. We have shown 
that the statistical behavior of the trapped trajectories 
is completely different from that of the free trajectories. 
The trapped trajectories have a quasi-coherent behavior. 
Their average displacement, dispersion and probability 
distribution function saturate. A very strong anomalous 
clump effect characterizes neighboring trapped trajecto- 
ries. Their clump life time is very large compared to the 
time of flight. This shows that these trajectories form 
structures similar with fluid vortices. The statistical pa- 
rameters of these structures (size, build-up time, disper- 
sion) are determined. The trajectories contained in such 
structures do not contribute to the large time diffusion 
coefficient. The later is determined by the free trajec- 
tories which have a continuously growing average dis- 
placement and dispersion. The probability distribution 
function for both types of trajectories are non-Gaussian. 

The time dependent diffusion coefficient is determined 
as a functional of the Eulerian correlation of the stochas- 
tic potential. 

The general conclusion of this work is that the exis- 
tence of an invariant in the evolution equation (the poten- 
tial) determines long-time correlations (memory effects) 
and coherence (trajectory structures) in the stochastic 
motion. 
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Figure 1 

Exemples of paths of average trajectories in subensembles (S2) obtained from Eq. (32) for several values of 4>ij and 

for the values of (fP and u mentioned on the figure. 
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Figure 2 

Exemples of paths of average trajectories in subensemblcs (S) obtained from Eq. (34) (continuous lines). The results 
obtained with the decorrelation trajectory method are also ploted for comparision (dashed lines). 
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Figure 3 

Time evolution of the average trajectory and of the dispersion [Eq. (35)] in subensembles (S) with = (a) and 

for 4>° = 1 (b). 
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Figure 4 

Contour plot of the pdf [Eq. (36)] of the trajectories in a subensemble (S) for {a.) (fP = 0, u = 1, t = ..Tfi and (b) 
0° = l,u = l,f = ..Tfi (at saturation). The path of the average trajectory is also represented (red line) on the 

interval [0,t] . 
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Figure 5 

The second moments of the relative displacement (^J(t))^ (continuous blue line) and (<5|(t))^ (continuous black 
line) compared to (^^(i))^ (dashed blue line) and (2^2 (0) 5 (dashed black line) for two subensembles (S): one with 
0° = 0, u = 1 (Fig. 5a) and the other with (f>" — 1, u — 1 (Fig. 5b). 
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Figure 6 

Contour plot of the pdf of the relative motion in a subensemble (S) for (a) cfP ^ 0, u = 1, t = ..Tfi and (b) cj/' — 1, 

u ^ 1, t = ..Tfi (at saturation). 
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Figure 7 

The time dependent diffusion coefficienr D{t) obtained from Eq. (38) with the nested subensemble method 
(continuous hne) compared with the result of the decorrelation trajectory method (dashed hne). 
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